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It has recently been suggested that the competition for a finite pool of microRNAs (miRNA) gives rise to effective 
interactions among their common targets (competing endogenous RNAs or ceRNAs) that could prove to be crucial 
for post-transcriptional regulation (PTR). We have studied a minimal model of PTR where the emergence and the 
nature of such interactions can be characterized in detail at steady state. Sensitivity analysis shows that binding 
free energies and repression mechanisms are the key ingredients for the cross-talk between ceRNAs to arise. 
Interactions emerge in specific ranges of repression values, can be symmetrical (one ceRNA influences another 
and vice-versa) or asymmetrical (one ceRNA influences another but not the reverse) and may be highly selective, 
while possibly limited by noise. In addition, we show that non-trivial coiTelations among ceRNAs can emerge in 
experimental readouts due to transcriptional fluctuations even in absence of miRNA-mediated cross-talk. 
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I. INTRODUCTION 

MicroRNAs (miRNAs) are 21-23 nucleotides (nt) long, 
endogeous, non-coding RNA molecules, that perform post- 
transcriptional regulation by specifically binding target mes- 
senger RNAs (mRNAs), typically leading to a reduction in 
the levels of the corresponding proteins |[T]-(3|. They are tran- 
scribed from independent miRNA genes or from introns of 
protein-coding transcripts. After being processed into ma- 
turity, a miRNA is loaded onto a specialized class of pro- 
teins to form the RNA-induced silencing complex (RISC), 
which specifically binds miRNA response elements (MREs) 
located in target mRNAs (usually in their 3' UnTraslated Re- 
gion or 3'-UTR) through a base-pairing recognition mecha- 
nism which requires at least 6-nt complementarity. The whole 
process, known as RNA interference (RNAi), results in gene 
silencing through translation inhibition and mRNA destabi- 
lization EHJ. 

Each mRNA can typically interact with several miRNAs, 
and each miRNA can target many different mRNAs. Within 
the complex network of potential interactions that ensues, 
miRNAs have long been thought to function mainly as fine- 
tuners for regulation by weakly dampening the protein out- 
put HI 121. This view is supported by the fact that the sta- 
tistically over-represented motifs (feed-forward or feed-back 
loops) that have been identified in the known miRNA-mRNA 
interaction network are indeed capable of buffering the noise 
level in the output layer (proteins) f6^-^\. More recently the 
attention has been directed to system-level effects. In par- 
ticular, it has been realized that miRNA-based regulation is 
strongly affected by global properties like the total concen- 
tration of available targets (a feature known as dilution effect 
[[101). In addition, it is now known that pseudo-genes and 
other long non-coding RNAs also possess MREs and can be 
bound by the RISC. This implies that, besides mRNAs, non- 
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coding RNAs sharing identical MREs compete for common 
miRNAs ifTTllTSll . Thus, miRNAs appear to mediate the cross- 
talk between a broad class of competing endogenous RNAs 
(ceRNAs), leading to a large-scale network of interactions 
across the transcriptome ifOlIT?! . Recent studies have shown 
that such interactions play a central role in many biological 
contexts, from muscle differentiation 1 12] to cancer lfT5l[T6l . 

Despite the body of experimental evidence, a clear quan- 
titative understanding of miRNA-mediated regulation is still 
lacking. To address this issue, we formulate a minimal model 
of post-transcriptional regulation and analyze its steady state, 
aiming at quantifying the intensity of the interactions aris- 
ing from competition through an analysis of the sensitivity 
to changes in the ceRNA transcription rates. We show that 
binding free energies and repression mechanisms are the key 
ingredients for the cross-talk between ceRNAs to arise. The 
emergent interactions can be symmetrical (one ceRNA influ- 
ences another and vice-versa) or asymmetrical (one ceRNA 
influences another but not the reverse) and may be highly se- 
lective, although possibly hampered by noise. 

We furthermore argue that the identification of cross-talk 
from gene expression data can be hindered by the fact that sta- 
tistically significant correlations among ceRNAs can emerge 
in the experimental readouts simply due to transcriptional 
fluctuations. 



II. RESULTS 

A. One miRNA species, N ceRNA species 

We start by considering a highly simplified system formed 
by different species of ceRNA molecules, labeled m, (/ = 
1 , . . . , A^), targeted by a single miRNA species labeled //. Each 
OT, can reversibly bind /i in complexes labeled c,. Conform- 
ing to the experimental evidence | IT] [TSji . we assume that 
translational repression is fast and precedes mRNA destabi- 
lization, implying that complexes cannot be translated and re- 
pression of translation simply occurs by sequestration of free 
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FIG. 1 Schematic representation of the considered processes for a 
system N = 2 ceRNA species and one miRNA 



ceRNAs. Taking a constant translation rate, the levels of un- 
bound ceRNA can be used as a direct proxy for protein con- 
centrations. The allowed processes with their respective rates 
are then as follows (see Fig. [T]i: 



d) ^ nii d)^ fi 11 + mi ^ ci 

dt 6 k- 







Ci Id 



Note that complexes are assumed to be degraded either 
through a catalytic channel (which gives back the miRNA to 
the cytosol) or through a stoichiometric channel (where both 
molecules are degraded with the complex). This choice serves 
the purpose of keeping the model as general as possible. The 
exact mechanism of target repression is still a matter of de- 
bate, and in the past years several mechanisms have been 
reported II2J3- Generically, miRNAs incorporated into the 
RISC do not seem to decay with their target, thus becoming 
again available for a new round of target RNA silencing. Nev- 
ertheless complexes may enter in specific cellular structure- 
like P-bodies thus resulting in an effective stoichiometric se- 
questration of both the miRNA and the target. Clearly, this 
setup represents a coarse-graining of the real biological pro- 
cesses, which typically requires multiple catalyzed elemen- 
tary steps (e.g. in the formation of the RISC). However, such 
details may be disregarded if, in presence of many different 
targets, the only rate-limiting factor is the miRNA concentra- 
tion. This will indeed be our main assumption (together with 
the fact that m,'s can only interact through fj). 

Denoting the concentration of species x by [x], we can write 



the mass-action kinetic rate equations for the above system as 



— [m,] = -di[mi] + bi - klUAim] + ^/ [q] 

= -m +p-Yj Kmm] + Yj^k + (1) 



— [q] = -(cr; + ki + Ki)[ci] + k^[M][mi] 

In turn, they lead to the steady state equations 

bi+kr[ci] 
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where we defined m* - bi/di, c* - bi/{cri + Ki) and 
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with (pi = kj^ I {(Ti + Ki). Note that m* and c* represent the max- 
imum achievable concentrations of free ceRNAs and com- 
plexes, respectively. A simple calculation shows that the bind- 
ing free energy of the complex (in units of kgT) is given by 



AG,- = - log 



k_t\mm 

kj{Ci] 



•log 
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(4) 



This clarifies the physical meaning of 0,-: if (pi » 1 then com- 
plex formation is close to equilibrium; if 0, <s; 1 instead, 
the process is unbalanced towards association. On the other 
hand, one sees that the quantity yUo,, defined in ([3]) (which 
only depends on the kinetic parameters of ceRNA m,) gives 
the miRNA level for which the concentrations of free ceR- 
NAs and of complexes equal half of their theoretical maxima. 
Therefore, in practice, we have the following situation: 

a. If [n] <K juo,/ (or [//]/;Uo,; = 0{e)), then [m,] ^ m* 
and [c,] <K c*: here the levels of free ceRNAs are 
largest, while complexes are roughly absent; sponta- 
neous degradation is the dominant channel of ceRNA 
decay. 

b. If [//] ^ jjoj (or |l - [fJ.]/noj\ = <9(e)), then [m,] ^ 
TO*/2 and c,- ^ c*/2: here free ceRNA concentration 
is roughly half the theoretical maximum; spontaneous 
ceRNAs decay and miRNA-mediated degradation have 
similar weight. 

c. If [ju] » yuo,, (or fioj/lp-] - 0(e)), then [m,] m* and 
[ci] ^ c*: here the levels of complexes are largest, while 
free ceRNAs are roughly absent; miRNA-mediated 
degradation is the prevailing channel of ceRNA decay. 

It is reasonable to expect that in cases a. and c. the steady 
state level of free ceRNA [m,] will only be weakly sensible to 
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FIG. 2 Top: steady state concentrations in a system withi N = 2 ceR- 
NAs, obtained by fixing all parameters but the transcription rate bi 
of ceRNA 1 . The shaded area highlights the difference between the 
steady state level of ceRNA2 with or without (i.e. for b\ = 0) its com- 
petitor mi: [nio] is sensible to variations of bi (via the change of the 
free miRNA concentration which is correspondingly dropping) 
only in an intermediate, narrow interval. Bottom: fractions of free 
molecules, namely (fii = [nJi]/([mi] + [ci]), (pi = [m2]/([m2] + [C2]), 
0fi = [/"]/([/'] + [t'l] + [C2]), versus bi. The dynamical range of the 
cross-talk interaction between the two ceRNAs corresponds to the 
window where free and bound molecules have similar concentration, 
i.e. to the <S regime. Parameters for this case (in their respective 

1, 



units) are as follows: ^2 = 10, rfi = A 
k: = kj = 10~^, a\ = (Tj = 10, ki = k2 = 



1, kl = 10\ ki, 



1. 



(small) variations in [//]. We shall call these regimes 'Free' 
and 'Bound' (iFand S for brevity), respectively. In case b., 
instead, the microRNA concentration lies in the dynamical 
range of F,-, so that [m,] will respond to (small) variations in 
{ji}. We will the call this regime 'Susceptible' {S for brevity). 
The outlook is that once the kinetic parameters of the ceRNAs 
are given, the are given as well, and the miRNA level 
suffices to know whether a ceRNA is in the ;F, >S or S regime. 
Note that analysis of the RNAi enzyme complex has shown 
that its kinetics can vary substantially across different targets 
and that it is strongly affected by the degree of complemen- 
tarity IIT9I. This means that, in principle, different targets may 
have very different Hofs and may thus be located in distinct 
regimes of regulation at fixed [//]. 

To illustrate the emergent interactions between ceRNAs, we 
plot in Figure |2] the steady-state levels of mi, m2 and yu as a 
function ofbi in a system with = 2 in which all other kinetic 
parameters are fixed. One clearly sees that, even within this 
basic model, a change in the transcription rate of a ceRNA can 
affect the steady state concentration of a different ceRNA. Our 
goal here is to characterize such a 'cross-talk' quantitatively 
by computing the 'susceptibilities' I, . . .,N) 



Xi. 
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which for / ^ j measure of the magnitude of the interaction 



between ceRNA m, and ceRNA nij. Using Q we get 

d[fi] dbj 



di 



(6) 



The cross-susceptibility (the last term in (j6])) can thus be seen 
as the product of two factors: the response of the miRNA level 
to perturbations of the transcription rate of ceRNA mj, and the 
response of the level of ceRNA m, ([m,] = m*Fi) to perturba- 
tions of the miRNA level. For the latter we get 



-l/fiQj ier 
-l/(Vo.,) ieS (7) 



Note that, expectedly, this function is negative and is largest 
for [fj] ^ //o,i, i.e. when the ceRNA is in the iS-regime. In 
order to compute d[^]ldbj, we need an explicit expression for 
[yu], which should be obtained from the steady state condition 
(|2]l. We re- write this as 



6 + Y^b■,z■,Fi{[^i\) 



(8) 



where z, = o-il{^iaj{cri + /c,)]. Equation ([8| teUs us that if 
cr, = for all / (i.e. if complex decay is purely catalytic) then 
[ju] -PIS and no cross-talk is achievable since [/i] is indepen- 
dent of bj. If however cr, + for some i, then other solutions 
are possible. In particular, (|8]l is an algebraic equation of or- 
der -I- 1 at most, an approximate solution of which can be 
obtained under the assumption that ceRNAs can be separated 
in the different regimes defined above. Using the fact that, up 
to next-to-leading order in e «: 1, 



Fm) - 



5-(M-/^o,,)/(4jUo,,) i^S (9) 
,y"o,i/[/^] i e S 



and neglecting (when necessary) terms of order 2 or higher, 
one finds that (see Supporting Text for a detailed derivation) 



P - ^iM'i - I 'Lies biWi 



with Wi = ZifiQ.i = o-ilicTi + Ki). One now sees that 



(10) 
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W, + M)/(W j^S (11) 
1 jeS 



where we have defined the shorthand 



Xtifi — 



dp 



ieT 



(12) 



Given a shift in the level of a ceRNA, the response of the 
miRNA is always negative (since an increase in bj causes an 
increase in the level of complexes cj). Also, if nij e T then 
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[fij/fioj <s: 1, i.e. the level of microRNA is roughly insensitive 
to small changes of the production rate of ceRNAs in the Free 
regime. 

Finally, combining (|7jl and ( 11 1 we obtain for Xij 



Xij = 



Fim)Sij 



(13) 



where Wr(,)_r(;) is a coefficient that depends only on the 
regimes R{i) and Ri j) to which / and j belong. In other words, 
R(i),R(j) e {T,.S,S) and the 3 x 3 matrix W is given by 



W 



^ y"0; [jU] I^Qj 



M 

4^ 



0(6) 



C(e) 
0(1) 
C)(e) 



0(e)^ 

0(1) 

0(e). 



(14) 

Three important observations can now be made about Xij- 
First, since both terms (|7]i and (111 are negative, the cross- 
talk between ceRNAs tends to correlate their levels. Second, 
the matrix W is not symmetric, as might perhaps have been 
expected. Finally, for / j all of the elements of W are of 
order e or smaller except Ws,s and Ws,s, which are of or- 
der 1 . This implies that, in this scenario, two types of effec- 
tive interactions arise: the first one encodes the response of a 
ceRNA in the ^S-regime to a perturbation of another ceRNA 
in the >S-regime, and it is symmetric; the second one encodes 
the response of a ceRNA in the »S-regime to a perturbation of 
a ceRNA in the S-regime, and it is not symmetric (i.e. per- 
turbing the 'susceptible' ceRNA the 'bound' one will not re- 
spond). 

In words, the scenario described here corresponds to a lin- 
ear response theory in which a change in the transcription rate 
of a ceRNA (i.e. bj bj + 5b j) induces a shift in [nij] (i.e. 
[nij] [nij] + Xjj^^i with xa > 0) and a shift in the level 
of miRNA (i.e. [ju] [/j] + Xtij^^j with Xfij < 0). In turn, 
this affects [m,] (i.e. [m,] [m,] +XijSbj with;^',^ > 0). So 
for instance if dbj > then [mj] increases, [fi] decreases and 
[m,] increases. The intensity of the interaction depends on 
the regimes to which the ceRNAs belong. In essence, high 
cross-susceptibility is achieved (i) (symmetrically) between 
ceRNAs in the >S-regimes, and (ii) (asymmetrically) from a 
ceRNA in the S-regime to a ceRNA in the .S-regime. The 
quantities fiQj cc l/kt can be seen to induce a hierarchy of in- 
teractions: ceRNAs in the S-regime (higher binding affinity) 
can unidirectionally affect ceRNAs in the ^S-regime, which in 
turn may influence other ceRNAs in the 5-regime. On the 
other hand, ceRNAs in the f regime (lower binding affin- 
ity) interact weakly with the rest of the system and fluctua- 
tions in their transcription rates do not propagate to other ceR- 
NAs. This defines a selective, possibly asymmetric channel of 
communication which links ceRNAs targeted by a common 
miRNA. It is important to remark that cross-talk appears only 
when the miRNA level is in a specific range, implying that the 
structure of the emergent interaction network is flexible and 
dynamical: the set of ceRNA species that interact may change 
upon varying [jj] . 




FIG. 3 Susceptibilities in a system of A' = 4 ceRNAs, as a func- 
tion of miRNA transcription rate y6 (all other parameters being fixed). 
In this example ceRNAs are cast in two groups: group A, formed by 
ceRNAs mi and nij, and group B, formed by ceRNAs m2 and mn. 
CeRNAs belonging to the same group share identical kinetic param- 
eters. In particular, /jq.i = y"o.3 ^ y"o.2 = Po,4- While for yS smaller 
than about 500 no cross-talk is observed, as /3 increases a symmetric 
interaction between ceRNAs in group A (of magnitude comparable 
to the self-susceptibilities) appears: X3i = Xa - Xa = Xn- As p in- 
creases further, this interaction is switched off, and ceRNAs in group 
B begin to cross-talk instead: = Xu - X22 = In this region, 
a change of transcription of a ceRNA in group A can affect the level 
of ceRNAs in group B, but not viceversa (asymmetric cross-talk): 
X2i = X23 = X4i = X43 » Xi2 = X32 = Xi4 = X34- Finally, for 
sufficiently large /3 all no cross-talk takes place. 



The emergence of selectivity and directionality as features 
of the cross-talk can be seen in a concrete case in Figure |3] 
where we plot the susceptibilities for a system of = 4 ceR- 
NAs. One sees that different interactions are switched on in 
different ranges of values for the miRNA transcription rate, 
leading to a gradual modification of the structure of the inter- 
action network as /3 changes. Notice that heterogeneity in the 
quantities /zo,i leads to interaction asymmetry. A schematic 
summary of the cross-talk in this system is given in Figure |4] 



We notice that the intensity of the cross-talk described by 
( 1 3 1 is modulated by the factor 



biWi 



4M "^MS + ZkesbkWk 



(15) 



where we used [/i] ^ yUo,* for k e S and neglected the con- 



tribution to III] due to ceRNAs in the ^-regime (see (lOl), 
which is of order e. Again, we see that an effective interac- 
tion requires some degree of stoichiometric degradation: if 
Wj = (i.e. if the complex decays in a purely catalytic man- 
ner) the cross-susceptibility vanishes. In addition, ( 15 1 sug- 
gests that the magnitude of the interaction is weakened only 
by ceRNAs lying in the >S-regime, so that, even in presence of 
a large number of interacting ceRNA species, cross-talk can 
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FIG. 4 Schematic representation of the patterns of interactions arising in a system of A' = 4 ceRNAs at different miRNA levels (increasing 
from A to D). A) All ceRNAs are in the T regime ([ju] <c /Jo.; Vi) and there is no interaction between them. B) ceRNAs 1 and 3 are in 
the S regime ([ju] ^ yUoj - yUo.s) and a symmetrical interaction between them is switched on. C) ceRNAs 1 and 3 are now in the S regime 
([ju] » ;Uo,i - /Jo.s) while ceRNAs 2 and 4 are in the S regime ([/j] ^ /Jo,2 - y"o,4): the resulting interactions are symmetric between 2 and 4 and 
asymmetric of ceRNAs 1 and 3 on ceRNAs 2 and 4. D) All ceRNAs are in the S regime ([/i] » //q,, V/) and no cross-talk occurs. 



be large if the overall population of ceRNAs in the >S-regime 
is restricted. On the other hand, many factors, such as the 
overall population of ceRNAs in the S-regime, affect the rate 
of miRNA transcription required in order for a given ceRNA 
to be susceptible. 

A rough approximate expression for the range AyS of val- 
ues of the miRNA transcription rate where ceRNA m,- is most 
responsive to the miRNA can be derived considering as 'sus- 
ceptible' a window of values such that /io.,72 < [/i] < 
3/io,i72. If so, then 



kWk 



(16) 



where we used ^ fio j. One sees that mainly depends on 
the transcription rate of the overall population of ceRNAs in 
the susceptible regime and on the degree of stoichiometricity 
of degradation. 

Finally, we observe that an analogous cross-talk scenario 
emerges for a system in which M miRNA species share the 
same target RNA: the level of a miRNA species may be highly 
susceptible to a change in the transctiption rate of a different 
miRNA when the level of the target RNA lies in a specific 



window (see Supporting Text for details). 



B. A' ceRNA species, M miRNA species 

Let us now consider the general case of a system formed by 
M miRNA species /Ua (a - 1, . . . , M) and ceRNA species 
m, (/ = 1, . . . , A^), defined by the rates 



bi 

^ rrii 



Pa Ka 

S„ k- 







and for which the following steady-equations hold: 



[nij] 
[Ma] = 



Pa + ILiiKa + Kia)[Cia] 



[Cia] = 



6a + 

klMa'W.fni] 

^ ia + ^ia 



(17) 
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Again, we shall focus on computing ceRNA sensitivities to 
perturbations of transcription rates of other ceRNAs. Ne- 
glecting higher order interactions involving two or more miR- 
NAs (which is justified if the miRNA-ceRNA network is suf- 
ficiently sparse) we have 



Xij 



d[mi\ 



d[mi] d[fia] 
^ d[^ia] dbj 



(18) 



The NM+N + M steady state equations ( 17 i can be reduced to 
N + M coupled equations for the N + M unknown m,) by 
eliminating the complexes. After some straightfoward alge- 
bra, we get [m,] = m*F,({[/i„])) and Ha = jU*^'a({['«i])) with 



where 



[mj] 

mo.ia 



di 



6a 



k:„, + Ki, 



CTla 



(19) 

(20) 
(21) 



In turn, the levels of free miRNA and ceRNA are described by 



[nii] 



5„ + ^ biZiaFi 
i 



^bi 



(22) 
(23) 



withZi^ = 0-ia/[flQjj0-ia+Kia)] and = (o-ia+KiJKmo^iaO-ia)- 

The quantity F, can be re-written as 



Fi^ 



1 



1 



Z<"^ 1 + [ju„]//io,« 



(24) 



where Z^"' = 1 + 2y^„[//y]/yUo,,y and /io,,a = f^OM^] 



(a) 



Note 



that sum in zf"' includes all miRNA species except for 



Equation (24i tells us that, in presence of many miRNA 
species, we may account for the effect of species on ceRNA 
m, by just re-scaling f , by z|"^ and shifting the reference level 
by zf"'. A simple interpretation of the above expressions 
can be gained by introducing an effective decay rate c/!"^ - 
diZ'f^ and noting that 



r 1 ^' 



Ao./ff = ^(1 



.) 



(25) 



One immediately recognizes the same form of the steady state 
equation v\ for the case M - \, and sees that Z?"* ultimately 
plays the role of a factor accelerating the effective turnover 

By analogy with the case M - 1, we will say that a ceRNA 
is free with respect to miRNA yu„ (and write / 6 Tia)) if 
[//q.] fio,ia\ it will be 'susceptible' with respect to ju„ (or 
; 6 S{a)) if \jia\ - fioja', it will be 'bound' with respect to 



fia (or / e !B{a)) if » /5o,m- Note that being bound with 
respect to a miRNA species is sufficient for a ceRNA to be 
translationally repressed. For consistency, a ceRNA can only 
be bound with respect to one miRNA species (in that case it 
will be free with respect to all other miRNAs). Separating the 
different regimes we have 



f[Zf^]-i(l-K]//io,,J 



/ e Tia) 



\ {ifWl - 0" - ma)mQM)\ i e S(a) 
In turn, for the levels of free miRNAs we obtain 



(26) 



\Ua] - 



^a + TjieTia) biZia + | TiieS(a) biZu 



(27) 



where Wia = (Tia I {(Tia + Kia) and Zia = 0-ial\fia,ia{cria + Kia)]- 

One may now compute the different terms of the susceptibili- 



ties. For the quantity (see ( 18 1) we finally get the M > 1 
analog of (pjj), i.e. 



Xi, 



1 

di 



Fidii + 



a;R(i),R(j) 



(28) 



where the matrices Wa are given by 



Wa = 



fJ^OJafJ^OJa 
[Ma]^ 



flO,ia fioja fioja 
iMa] fioja + \Ma] [fia] 
fioja 
jioja 



fio,iaP-OJa f^OJa 4/iojQ' 
. fioja flOJa iioja + [fJ^a] 



fioja [fJ^a] fioja 



[f^a] 



(29) 



In this case, the intensity of the cross-talk described by ( 28 1 is 
modulated by the factor 



biWjaXaa 



biWj 



AZf'^fia] Z^^'iAljjM. + HkeSia) bkWka) 



(30) 



We therefore conclude that miRNA yu„ gives a relevant contri- 
butions to the overall susceptibility Xij if either: 

(i) / e Sia); 

(ii) j e S{a) or j e S(a) (in the latter case, fi^ is the main 
repressor of my); 

(iii) z'"* ^ 1, i.e. ceRNA m,- has few repressors besides 

(iv) few ceRNA species belong to S(a), so that dilution is 
limited. 

In summary, the effect of 'background' miRNAs which do not 
mediate interactions is an increase in the effective rate of de- 
cay, and consequently a shift in the susceptibility threshold. 
On the other hand, the effect of background ceRNAs is a di- 
lution of the cross-talk among ceRNAs, as seen in the case 
M - \. An illustrative example of the interactions arising 
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FIG. 5 Schematic representation of a system of A' = 4 ceRNAs 
species and M = 3 miRNA species. Continuous blue arrows link 
miRNA a to ceRNA i if i 6 Bia); dashed red arrows are found 
if / 6 <S(q'); dotted green arrows are found if / 6 T'(a). In this 
case, ceRNAs 1 and 3 are both in S{A) and S{B), ceRNAs 3 and 
4 are in S(C), and ceRNA is 1 in S(A). This situation results in 
the following interactions: symmetric cross-talk between ceRNAs 1 
and 3, mediated by miRNAs A and B; symmetric cross-talk between 
ceRNAs 3 and 4, mediated by miRNA C; asymmetric cross-talk from 
ceRNA 2 to ceRNA 1, mediated by miRNA B; asymmetric cross-talk 
from ceRNA 2 to ceRNA 3, mediated by miRNA B. 



among ceRNAs in a system with N - A and M = 3 is shown 
in Figure |5] 

To conclude, we notice that network topology can play an 
important role as interaction enhancer. For instance (see Sup- 
porting Text), cross-talk can take place among ceRNAs in the 
Free regime (in spite of the small ;^', y;„) provided they are com- 
monly targeted by a large number of miRNA species. In other 
terms: interactions between ceRNAs can be mediated by a 
large number of miRNA species which individually would 
only weakly dampen ceRNA levels. However, in order to 
achieve efficient cross-talk strong correlations in the network 
connectivity are needed, so that highly clustered networks can 
allow for much stronger cross-talk than random graphs. 



C. Steady state fluctuations 

Genetic circuits that regulate cellular functions are subject 
to stochastic fluctuations, specifically in the levels of the dif- 
ferent molecular species that interact ll20l[2T]| . Noise, far from 
being just a nuisance, plays an essential role in cellular activi- 
ties, for example by enabling coordination of gene expression 
across large regulons, or by allowing for probabilistic differ- 
entiation of otherwise identical cells ["221. On the other side, 
a noisy gene expression is potentially harmuful in many sit- 
uations: in developmental circuits, for example, it can lead 
either to arrested development, aberrant positional expression 
of tissue specific genes or over-representation of specific cell 
types |23|. If there is only a relatively narrow protein level 
which is optimal, some sort of tuning must act to prevent fluc- 
tuations outside the functional range. Cross-talk of the type 
discussed so far may either result in an amplification of up- 
stream fluctuations or represent an efficient noise buffering 
mechanism. To analyze this issue in some detail, we focus 
on the role of transcriptional noise, the primary cause of vari- 



ability in gene expression among cells in isogenic populations 
||24j. If one assumes that extrinsic transcriptional noise is the 
dominant source of stochasticity and neglects molecular noise 
entirely, it is possible to estimate concentration fluctuations 
in the ceRNA-miRNA networks at steady state, obtaining ex- 
pressions valid in the linear response regime. 

Let us consider for simplicity a system of N - 2 ceRNA 
species and M - \ miRNA species, and let P(r) denote a dis- 
tribution of transcription rates (where r = {bi,b2,/3}), such 
that an ensemble of systems at steady state can be constructed 
by sampling a vector r from P(r) for each system in the en- 
semble. For P(r) one may for simplicity take a Gaussian, i.e. 



P(r) cc exp 



-^(r-r)^2-Hr-F) 



(31) 



where r - {bi,b2,/3] is the mean and S is the correlation ma- 
trix of inputs. Clearly, a distribution of transcription rates in- 
duces a distribution of steady state concentrations. The latter 
is what we aim at characterizing. 

If variability in transcription rates is sufficiently small, we 
can expand the steady state levels ( - {[mi], [022], [ju]} (note 
that { = ({r)) around r (small noise expansion), obtaining 



£i - + 



k 



Xik 



(32) 



where = /',(r). In this approximation: 
Prob(f = x) = J" P(r)S[e(r) - x]dr = 



— Nexp 



(33) 



where X - (;f-"')^2"'^"', x being the matrix of susceptibili- 
ties defined in ( 32 1. The joint probability distribution and the 
susceptibility matrix can then be used to characterize steady 
state fluctuations and correlations, e.g. 



0-] = Hi 



YjXijXik^jk 

j.k 



(34) 



For uncorrelated transcription rates the covariance matrix E 
is diagonal and (|34j) reduces to crj - YjkX'fk'^kk' expectedly, 
each term positively contributes to increase the noise. As 
shown in Figure |6] fluctuations can become very large in the 
susceptible regime as the system is strongly coupled, possibly 
limiting the efficiency of signaling in ceRNA network . 

In presence of coiTelations at the transcriptional level, how- 
ever, the signs of off-diagonal terms become crucial. Re- 
calling that, generically, Xi,ii < 0^ Xij > 0' ™d XfiJ < 0' 
one sees that anti-correlated ceRNA transcriptions and corre- 
lated miRNA-ceRNAs transcriptions may lead to a reduction 
of fluctuations with respect to the uncoiTelated case, as shown 
again in Figure|6] On the other hand, negative miRNA-ceRNA 
coiTelations and positive ceRNA-ceRNA coiTelations strongly 
amplify fluctuations. In other terms, miRNA-mediated cross- 
talk coupled with correlation of transcriptional inputs may 
represent a powerful noise processing mechanism. 
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FIG. 6 Level fluctuations induced by transcriptional noise in a sys- 
tem with N = 2 ceRNA species and M = I miRNA species. We 
took Gaussian distributions for the transcription rates /?, b] and b2, 
keeping the ratio between the average and the width fixed in each 
case. All parameters and distributions of rates are also kept fixed, 
except for the miRNA transcription rate distribution P(J}), which is 
parametrized by its average value /J. In this case: b[ = b2 = 10', 



lo^ kv 



0, a-: 



10, 



K\ = K2 = 0. Top: normalized fluctuations (ratio between the width of 
the fluctuations in the interacting and the non-interacting system, the 
latter corresponding to = = 0) for uncorrelated distributions 
of transcription rates. Center: ratio between the normalized fluctu- 
ations of the level of ceRNA I obtained in presence of correlations 
and in the uncorrelated case. Yellow line: maximal anti-correlation 
between b^ and bi (J^n = ~^)\ purple line: maximal correlation be- 
tween bi and /? (E13 = 1); blue line: maximal correlation between ^2 
and p (E23 = 1)). Bottom: average molecular levels. 



D. Detection of miRNA-mediated cross-talk from gene 
expression data 

A key issue of the ceRNA scenario concerns the detection 
of cross-talk in gene expression data, typically from coiTela- 
tions or related quantities. It is important to note that, within 
the theoretical framework we discuss, the presence of statis- 
tically significant correlations between ceRNAs is not neces- 
sarily a signature in this sense. Indeed, the Pearson correla- 
tion coefficient between ceRNAs, for independent transcrip- 
tion rates, reads 



P12 - 



TjkXikXikcrl 



(35) 



where (t\ is the variance of r^.. However \f xn - then 



Pi2 =A 



5[mi] <9[m2] 



(36) 



with A > a constant. Since both susceptibilities on the right- 
hand side are negative, a positive correlation between ceRNAs 
can emerge also in absence of miRNA-mediated cross-talk. 



More recently, information theoretical quantities have been 
employed as a means to detect miRNA-mediated cross-talk. 
In lfl4ll . for instance, the functional 

A/([mi], [m2]) = <[/([mi], [^i])]\,„^^ - I{[m,], (37) 

(I(x, y) denoting the mutual information of random variables x 
and y, {■ ■ ■ )z denoting the average with respect to the random 
variable z) has been proposed, with the rationale that if A/ > 
then the knowledge of [mi} increases the mutual dependence 
of [mi] and [/z], which can be interpreted as a signature of 
cross-talk between m\ and m^- We note, however, that under 
the steady state equations ([TTji free ceRNA levels depend only 
on the levels of the miRNA they interact with, so that the joint 
probability distribution of the levels of the various molecular 
species involved can be factorized, e.g. for = 2 and M - 1 



P([mi], M, M) = Pi([mi], bu])P2(N2], M) 



(38) 



This in turn implies that the three-species correlation func- 
tions can also be factorized, leading to A/ = independently 
of there being cross-talk or not. 

Quite importantly, however, in a typical experimental out- 
put (e.g. by microaiTay or deep sequencing analysis) it is hard 
to disentangle the contribution of free and bound ceRNAs and 
the experimental readouts give proxies for the quantities 

[mihp = [mi] + [a] , [fihp = [ju] + [c,] . (39) 

Based on Q, one has in particular 



[mi]xp = c- 



di 



l+2c*F,(M)/Mcy 



(40) 
(41) 



Note that if the lifetime of complexes (cr, -1- a-,) is shorter than 
that of ceRNAs d,, as is reasonable to expect [m,]^p, like 
[mi], is a decreasing function of [/i]. One sees that, in gen- 
eral, it is not possible to express the experimental ceRNA lev- 
els in terms of the miRNA levels only. Therefore the quan- 



tity ( 37 1 computed using the experimental readouts ( 39 1 can 
be different from zero. Again, however, this is not neces- 
sarily a signature of cross-talk. An argument is given in 
Supporting Text, where it is shown how non-zero values of 
AI{[mi]j:p, [fJ.]xp', [m2]xp) can be obtained even in absence of 
stoichiometric complex degradation (and hence of cross-talk 
at stationarity). 

In summary, more refined detection methods are likely to 
be needed in order to identify cross-talk among ceRNAs from 
gene expression data. 



III. DISCUSSION 

Recent experimental studies have suggested that the 
miRNA-mediated competition between ceRNAs could con- 
stitute an additional level of post-transcriptional regulation, 
playing important roles in many biological contexts. Trying to 
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achieve a clear quantitative understanding of the emergence of 
this effect has been the goal of this work. We have presented 
a minimal, rate equation-based model that is able to describe 
the cross-talk arising from competition at steady state through 
a systematic analytical characterization of the sensitivity to 
small changes in the transcription rates. To keep mathematical 
complexities to a minimum, we have adopted a coarse-grained 
view of the real biological process, even neglecting details of 
the miRNA-mediated regulation that could impact the emer- 
gence of cross-talk among ceRNAs. For instance, binding 
to the Argonaute/Ago protein (the catalytic component of the 
RISC) may represent a significant rate-limiting step |25|, and 
the competition for Ago has been shown to contribute to the 
emergence of ceRNA-ceRNA cross-talk |26|. 

The emerging scenario, valid in the linear response regime, 
is rather rich and complex. Interestingly, the competitive in- 
teractions can give rise to a rather selective communication 
channel: only ceRNAs in an intermediate, susceptible regime 
are responsive to miRNA perturbations and significantly con- 
tribute to diluting the strength of the interaction. Thus, 
even in case of a dense miRNA-ceRNA network, the result- 
ing ceRNA-ceRNA cross-talk pattern may be rather sparse. 
Moreover, interactions switch on only in specific ranges of 
miRNA concentrations, so that the structure of the emergent 
ceRNA-ceRNA network can adjust in response to variations 
in the miRNA levels. Perhaps unexpectedly, heterogeneity of 
kinetic parameters can give rise both to symmetric and asym- 
metric couplings. Furthermore, an analogous cross-talk sce- 
nario emerges between different miRNA species sharing the 
same target RNA. And, finally, the topology of the ceRNA- 
miRNA network may play an important role as strong correla- 
tions in connectivity in that network can enhance the ceRNA- 
ceRNA cross-talk. 

The above picture requires that miRNA-ceRNA complexes 
decay, at least partially, through a stoichiometric channel of 
degradation: for purely catalytic decay no cross-talk is possi- 
ble at stationarity, and perturbations of transcription rates only 
cause a transient response. Dynamical effects may neverthe- 
less play an important role on the time-scales of many cellular 
processes, and will be explored in a forthcoming work. 

In order to evaluate the robustness of the miRNA-mediated 
coupling, we have also performed a basic analysis of the im- 
pact of noise. Assuming extrinsic transcriptional noise as the 
dominant source of stochasticity, we estimated level fluctua- 
tions in the ceRNA-miRNA networks at steady state, again 
obtaining expressions valid in the linear response regime. It 
turns out that miRNA-mediated cross-talk, coupled with cor- 
related transcriptional inputs, represents a powerful noise pro- 
cessing mechanism that can lead to either noise reduction or 
amplification. It is interesting to observe that a circuit display- 
ing specific transcriptional correlations has been discussed in 
lfT2J . where a muscle-specific miRNA (miR-133b) embedded 
in a non-coding transcript (linc-MDl) has been identified. 
Clearly, linc-MDl 's transcript acts as a very efficient decoy 
for miR-133b. A theory for this case is worked out in the Sup- 
porting Text. 

It would be important to carry the analysis of the role of 
noise beyond the steps discussed here. Post-transcriptional 



regulation based on stoichiometric repression has been shown 
to cause large intrinsic fluctuations in intermediate regimes 
of repression |27|, effectively posing a limit to the possibility 
of having an efficient quantitative signaling between ceRNAs. 
Our analysis suggests on the other hand that cross-talk medi- 
ated by large number of miRNAs might be more robust. A 
more thorough mathematical/computational analysis, includ- 
ing molecular noise, may be able to shed light on this impor- 
tant aspect. 

We have finally shown that non-trivial correlations among 
ceRNAs can emerge in experimental readouts due to tran- 
scriptional fluctuations even in absence of miRNA-mediated 
cross-talk. 
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SUPPORTING TEXT 

A. Derivation of Equation pO) 

An approximate, explicit expression for the steady state 
miRNA level can be derived starting from (|9]). For simplic- 
ity, let us consider the case of N - 2 ceRNAs, in which (|8]l 
reduces to the equation 



[//] [6 + b,ZlF^{[^i]) + b2Z2F2{[^l])] = p 



(42) 



We can work out its solutions explicitly depending on the 
regimes to which the ceRNAs belong by inserting (j9]l into ( 42 1 
and keeping only linear terms in [//]. One finds the following 
results: 



(a) mi,m2 : 



(b) nix, 1112 G S 



5 + Hi ^1^1 



(c) mi,m2 & S : 



(d) mi eS, m2eS : 

W ~- 

(e) mi eT, m2 g S : 

(f) mi e T, m2e S : 



P - biW\ - W2 

S + \b2Z2 



P - b2W2 

6 + biZi 



P - \b2w2 



(43) 



(44) 



(45) 



(46) 



(47) 



(48) 



5 + bizi + \b2Z2 
Extending to the general case of ceRNAs we conclude that 

P - Hies biWi - \ Y,ieS biWi 



(49) 



B. The dual system: one ceRNA species, M mlRNA 
species 

The dual system in which M miRNA species, labeled 
{a - 1, . . . , M), target the same RNA m (to avoid confusion 
we will keep referring to it as a ceRNA even though in this 
case it is not really competing, being the only target species). 



can be worked out in full analogy with the case discussed 
above. In particular, defining fi* = Pa/^a, we have that, at 
stationarity, the level of free miRNA species is given by 



iMa] = H*aF,([m]) , 



where 



Fa 



[m] + mo,a 



mo,a 



da 



(1 + If/a) 



(50) 



(51) 



with ifTa - (k^ + Ka)/cra- (Note that now rates carry the index 
of the corresponding miRNA involved.) The free ceRNA level 
on the other hand results from the algebraic equation 



[m] 



d + Y,Pa^aFa([m]) 



(52) 



where - (cTa + Ka)/(niQ aO-a)- As before, each niQ a can be 
interpreted as reference level for the target which can be used 
to separate different regimes for the miRNA species. Borrow- 
ing the terminology used in the previous case, we have 



Fai[m]) 



1 - [m]/mQa a eT 

\ - ([m] - mo_„)/(4mo,„) aeS 
mo,a/[m] Of e S 



(53) 



where [m] <s niQ a for a 'free' miRNA, [m] ^ m^a for a 'sus- 
ceptible' miRNA, and [m] » rtio a for a 'bound miRNA. In 
turn, the level of free ceRNA is given by 



[m] 



(54) 



where Ua - ^a '«o,a = (o"ff + Ha) la a, while for the susceptibil- 
ity we obtain 



Xay — 



d\P-a\ 

dPy 



1 

Sa 



Fa{[m\)Say + 



m, 

4[m] 



-W, 



R(a),R(y) 



(55) 



d[m] 



and the matrix W is the same as in with 



where Xm,m - 

[to], too.q. and toq ^ replacing respectively [yu], /io,i and yUoj- 

Therefore the cross-talk that is established between miR- 
NAs is, as before, selectively turned on only for species lying 
in particular regimes, defined by the free ceRNA level. In 
complete analogy to the dual system analyized in the main 
text, two types of effective interactions arise: the first one is 
symmetric encodes the response of a miRNA in the >S-regime 
to a perturbation of another miRNA in the >S-regime; the sec- 
ond one is asymmetric and encodes the response of a miRNA 
in the >S-regime to a perturbation of a miRNA in the S-regime. 
An example of a pattern of interactions between miRNAs is 
shown in Figure |7] Note that the intensity of the cross-talk 
is modulated by the factor (y and increases when the rate of 
catalytic degradation increases. If cr^^ —» for all a implies 
'«{).« — > in this case, all miRNA species lie in the ;F-regime 
and no cross- talk is possible at steady state. 
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^1 miRNA 
Q ceRNA 




H H H H H 




FIG. 7 Schematic representation of a system of one target RNA and 
M = 6 miRNAs species. In this case miRNA 3 is Bound, miRNAs 
2 and 5 are Susceptible and the remaining are Free from the target 
RNA. Cross-talk interactions pattern is derived analogously to the 
dual case discussed in the main text: symmetrical cross-talk interac- 
tions emerge between miRNA 2 and 5 and asymmetrical interactions 
emerge from miRNA 3 to miRNAs 2 and 5. 



C. The role of topology 

Network topology can play an important role as a cross- 
talk enhancer. In specific, we will now argue that ceRNA- 
ceRNA interactions can be mediated by a large number 
of miRNA species which individually would only weakly 
dampen ceRNA levels. We consider a diluted network de- 
scribed by an adjacency matrix {Aja} such that 



A;^ — 



1 if ceRNA m,- is targeted by miRNA 
otherwise 



(56) 



making the following simplifying assumptions: (a) the net- 
work is kinetically homogeneous, i.e. rates are the same for all 
ceRNAs, so that jjQ^ia = jUq for each / and a and bi - b, d, - d 
for each (b) miRNA levels are uniform, i.e. - [//] for 
all a; (c) [fj] / fiQ = t <& 1 , so that all ceRNA species are in the 
;F-regime with respect to any miRNA (/ e 'Fia) V/, a). 
Consider a pair of ceRNAs m, and nij, targeted respec- 



tively by «/ 
n 



Yja^ia and rij - Yja-^ja miRNA species. 

In this case, the 



Yja-'^ia-Aja of which are in common 



ceRNA concentration reads m; - m*/(l+«,f) (withm* = b/d) 



and the cross-susceptibility ( 28 i turns out to be given by 



Xij.a 



lo 



if 



A- A ■ - 1 
otherwise 



(57) 



Consider now the particular case of a regular bipartite net- 
work with fixed ceRNA and miRNA connectivity so that 
n, = n for each ; and v„ = 2, A/^ = v for each a. Setting 



6 V 

Ka^K = — + 

zb I +tn 



(59) 



for all a we clearly see that now each miRNA species con- 
tributes equally to the overall susceptibility, i.e, Xij.a = Xo for 
all ( and j targeted by /Iq, with 



Xo 



1 



dK [1 +t(n - 1)][1 +tin- 1)]2 



(60) 



while the overall susceptibility is given by xu - nijXo- The 
contribution of a single miRNA to the overall susceptibilities 
will depend on the value of t. In particular, one easily sees 
that 



Xa 



'0{'-) forf<s:l/« 
:k~OC-) forf^l/« 

0{'-) forf»l/n 



dK 
1 

dKn 
1 

^ dKfiir' 



(61) 



Generalizing the Free, Susceptible and Bound regimes, one 
realizes that the case t l/n (resp. f ^ 1/n and f » 1/n) de- 
scribes a ceRNA that is 'globally free' (resp. 'globally suscep- 
tible' and 'globally bound') with respect to the overall miRNA 
population. We therefore conclude that Xij 

(i) increases with the number «y of miRNA species shared 
by the ceRNAs m, and my, 

(ii) decreases if the shared miRNAs have many other tar- 
gets; 

(iii) peaks when ceRNAs are 'globally susceptible' to the 
overall miRNA population, and it can be of the same or- 
der of magnitude as the self-susceptibility, i.e. 0(1 /d), 
when tiij ^ n. 

Perhaps most remarkably, the cross-talk can be effective even 
among ceRNAs that are in the Free regime with respect to in- 
dividual miRNAs, provided they are commonly targeted by 
a large number of miRNA species thus becoming 'globally 
susceptible'. However, in order to achieve efficient cross- 
talk strong correlations in the network connectivity are needed 
(large rtjj): highly clustered networks can allow for much 
stronger cross-talk than random graphs (see Figure[8|. 



where Kc - [S/(zb) 
the fact that zia - 



(o-+/f)//o 



Ztea(l + f«i) '] and z is defined by 
- z for each / and a. (The notation 



k e a indicates all ceRNAs interacting with miRNA /Zq,.) As 
expected, the dilution increases upon increasing the number 
of ceRNAs interacting with a given miRNA species /v„ (each 
of them add a positive term (1 -i- fnt)"' to thus making it 
larger) and upon increasing n, and rtj, since 



1 



Xi, 



(ni,nj » 1/r) 



(58) 



D. The mlRNA-decoy transcript 

Many miRNAs (possibly about 50% of the total H]) are 
hosted in non-coding genes whose transcript can incur a dual 
fate: after transcription, the precursors can either be processed 
into mature miRNAs through a series of steps involving pro- 
teins DROSHA and DICER, or they can reach the cytoplasm 
unprocessed in the form of long non-coding RNAs (IncR- 
NAs). The RNA sequence close to the sites corresponding 
to the miRNA presents a region with a sequence that is almost 
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(2) 






B 



FIG. 8 Two examples of different network structures with N = 2 
ceRNAs (blue circles) and M = 1 miRNAs (red squares). A) A 
highly correlated network structure where ceRNAs share almost all 
of their regulators (ni = ni = 5, nu = 4). B) A poorly correlated 
structure where ceRNAs share a small fraction of their regulators 
(«i = 112 = 4, ni2 = 1). Cross-talk will tipically be much stronger in 
A than in B. 




r(1-a, 



miRNA |j 



m (decoy) 



FIG. 9 Schematic representation of the model of a miRNA-decoy 
transcript. 



complementary to that of miRNA. These proximal strings al- 
low for the miRNA precursor (pri-miRNA) to take on the pe- 
culiar hairpin structure that is essential for the recognition by 
the processing proteins and thus for miRNA maturation fj)- 
It also follows, however, that the RNA sequence close to the 
miRNA necessarily contains a good potential binding site for 
the miRNA itself. When matured into IncRNAs, such tran- 
scripts are thus targeted by the miRNA and represent efficient 
'miRNA traps' or decoys, through which the population of 
miRNAs available for target repression can be regulated. The 
above miRNA-decoy mechanism can be modeled with follow- 
ing processes (see also Fig. |9]l: 



„ b 



q 



f(l-Q-) 

q ^ H 



(62) 



including transcription of the long non-coding RNA q at rate 
b, transport of q to the cytoplasm with processing into mature 
miRNA yu at rate (1 - a)r, and transport of q to the cytoplasm 
ar. The quantity 1 - a e [0, 1] thus gives the fraction of 
miRNA produced over the total number of transcribed RNAs. 




FIG. 10 Pearson correlation coefficient between the production rate 
of miRNA /j and of decoy m, for different values of the processing 
noise level (A = (r^^l[a( \ — a)} on the x axis) and of the transcription 
noise level (S = crl/b- on the y axis). High level of processing noise 
gives rise to negative correlations, while low level of processing noise 
and high level of transcriptional noise result in positive correlations. 



At stationarity, the miRNA and the IncRNA m are produced 
at constant rates according to 

m — ba (63) 
fi^b(l-a) (64) 

If noise affects both the transcription rate b and the processing 
efficiency a (taking again Gaussian distributions with means 
b and a and variances cr^ and cr^, respectively), the covariance 
between production rates is easily seen to be given by 

mfi- m ft - crl(a - a^) - o-^^b^ (65) 

Hence noise in b and a induces noise at the level of molecu- 
lar concentrations, yielding either positive or negative correla- 
tions between the steady state production rates of the miRNA 



jj. and of the decoy m as shown in FigurefTO] (Clearly, this con- 
clusion holds as long as the noise on a is sufficiently small, or 
A is not too close to 1.) These correlations, in turn, can result 
in a change of steady state fluctuations of other competing 
RNAs through the usual miRNA-mediated channels. In the 
case of muscle differentiation discussed in |12|, large levels 
of noise at the transcriptional or at the processing level could 
be exploited in order to increase cell variability and give lise 
to the differentiation program. Such a mechanism could be 
shared by other miRNA genes representing a widespread net- 
work motif. 
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E. On the significance of the conditional mutual 
information as a means to signal cross-tall< 

Consider a system (t,m,fi) of 2 ceRNAs (a target f and a 
modulator m and background targets) and one miRNA fj., 
subject to transcriptional fluctuations. Let us say that the ex- 
perimental readouts concern the quantities 



[m]_^p = [m] + [Cm] 
[t]xp = [t] + [c,] 
= [//] + [c,„] + [c,] 



Ixp 



(66) 
(67) 
(68) 



where [c,] and [c,„] represent the levels of miRNA-target 
and miRNA-modulator complexes, respectively. Suppose that 
both complexes decay catalytically, i.e. that the rates of stoi- 
chiometric complex degradation cr,„ = cr, = 0. In such condi- 
tions no cross-talk is possible at steady state. Furthermore, let 
us assume that the transcription rates b,, b,„, and /3 are drawn 
from a probability distribution Poibi,b,„,/3) such that 



Po(b„b,„,(^) = P{b„mbm-k) 



(69) 



with P an unspecified probability distribution with finite 
covariance (i.e., that the target and miRNA transcription 
rates are random variables while the modulator transcrip- 
tion rate is fixed at k). We want to show that, in 
this case, A/([f] ^p, [iJ.]xp', [>n]xp) > (with A/ defined in 
( 37 1) necessarily. This would imply that the condition 



^Klt]xp, [lJ]xp', [m]xp) > cannot be considered as a suffi- 
cient condition for cross-talk, since knowledge of [m] ^p can 
increase the mutual dependence between [filjcp and [t] ^p even 
in absence of cross-talk. 



To see this, note that the measured steady state levels are 
stochastic variables which depend on the transcription rates 
as 



[m]xp = fmd^) 
Mxp = Mb„P) 
[t]xp = Mb„/^) 



(70) 
(71) 
(72) 



(with f„,, ff, and f, unspecified functions). Now let us focus on 
(70 1 and (71 1. Given their monotonicity with respect to each 



of the variables on which they depend, they can be inverted: 



/3 = .CiMxp) 

b, = fp'mxp,ii) 



(73) 
(74) 



Hence it is possible to express [f] as a function of [m] and 
directly: {t\xp - h{{m'\xp,{lAxp)- In other terms, one finds 
a deterministic dependence of {t\ ^p on {m]xp. This implies that 
for each fixed [m]^^ the mutual information between [t]xp and 
{^i}xp diverges. As a consequence, their mutual information 
averaged over [m] ^p, \Ii{t]xp, [^Axp)) . diverges as weU. At 
the same time, however, the mutual information between {t\xp 
and [p]xp stays finite due to the noise on b, and yS. Hence 



imxp, Mxp)), , -mxp, Wxp) > 

(75) 

necessarily. 



